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Abstract 

We report the experimental observation of the photothermal effect. The measurements are 
performed by modulating the laser power absorbed by the mirrors of two high-finesse Fabry-Perot 
cavities. The results are very well described by a recently proposed theoretical model [M. Cerdonio, 
L. Conti, A. Heidmann and M. Pinard, Phys. Rev. D 63 (2001) 082003], confirming the correctness 
of such calculations. Our observations and quantitative characterization of the photothermal effect 
demonstrate its critical importance for interferometric displacement measurements towards the 
quantum limit, as those necessary for gravitational wave detection. 

PACS numbers: 78.20.Nv, 72.70.+m, 04.80.Nn 
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The study of thermal noise in optical components has a great interest both for the basic 
physics of the involved phenomena and for its implications in determining the fundamental 
limits of displacement measurements. In particular, the surface of the mirrors undergoes 
fluctuations which can depend on several dissipation mechanisms. Brownian noise, or vol- 
ume fluctuations, is generally related to homogeneous dissipation inside the mass of the 
mirrors, and is studied long since However, other dissipative mechanisms can in- 

duce surface fluctuations which can exceed the Brownian ones. Recently several authors 
examined the effects of temperature fluctuations which translate into surface fluctuations 
via the thermal expansion coefficient |3], f|, U; in particular they considered thermodynamic 
fluctuations of temperature of the mirrors and fluctuations caused by randomly absorbed 
photons (photothermal effect). Thermodynamic fluctuations of temperature are an intrinsic 
effect of the mirror material (substrate and reflecting coating) and are independent of the 
radiation power on the mirror. On the contrary, photothermal effect, beside the dependence 
on the thermal property of the mirror materials, depends on the power impinging on the 
mirror surface as well. Indeed, when light reflects off the surface of a mirror, a fraction of 
it is absorbed at the surface causing a local temperature change and a local deformation of 
the mirror surface. A fundamental limit to the position fluctuations due to photothermal 
effect is given by the standard quantum noise of the laser. 

Accurate studies of thermal effects are of crucial importance for the applications to grav- 
itational wave detectors. The sensitivity of long arm Michelson interferometers || depends 
critically from such effects, and also for massive detectors the most recently studied readout 
schemes rely on opto-mechanical techniques 0. In particular, Braginsky et al. || pointed 
out that the performance of advanced interferometric GW detectors, such as LIGOII [|J, 
might be limited by thermodynamic fluctuations of temperature and by photothermal effect. 
They derived an approximated expression for the spectral power density of photothermal 
noise which is valid for time scales smaller than the thermal relaxation time r c (adiabatic 
limit), i.e., for angular frequencies u larger than 



where k is the thermal conductivity of the mirror, p is the density, C is the specific heat 
capacity and r is the beam radius (assuming a Gaussian transverse profile), defined as the 
radius at which the power drops to 1/e of its central value. More recently the calculation 
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has been extended by Cerdonio et al. || to angular frequencies u smaller than the adiabatic 
limit, defined by Eq. ([I]). They predicted an expression for the power spectral density of the 
mirror surface position Sp T due to the photothermal noise, as 
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s PT (u) = 4 (i + °Y % K(n) , (2) 

where a is the Poisson ratio, a is the thermal expansion coefficient, S a \y S (uj) is the spectral 



power noise of the absorbed light [ 14 1 , and 

2 



1 f°° r°° r u 2 e" a2/2 
du I dv 



(3) 



n Jo J-oo (u 2 + v 2 )(u 2 + v 2 + iQ) 

is a shape function of the adimensional variable Q = uj/uj c . Assuming a flat S'absj as in the 
case of shot noise, this model predicts that for Q ^> 1 Spt is proportional to Q~ 2 , whereas 
for Q <C 1 it is almost constant. 

There exist some experimental evidence that the mirror heating changes the lenght of 
optical cavities. In Ref. || the authors describe the measurement of the output power of 
a CO2 laser at the turn-on and find a good agreement with simulations of the laser opti- 
cal cavity where the mirror surface is deformed by heating. In a different experiment, An 



et al. [|10j measure the transmission through a high-finesse Fabry-Perot when the laser is 
rapidly swept. They find an hysteretic behavior well explained by considering the mirror 
heating. Though their observations represent an indication of photothermal effect, their 
measurements are limited to the adiabatic limit defined by Eq. (jl|) and can hardly be com- 
pared with accurate models of the frequency response. So far the mentioned theoretical 
predictions have not yet been verified by experimental investigations, which are indeed the 
subject of the present work. 

In this paper we report for the first time on the experimental measurement of the fre- 
quency behavior of the photothermal effect: the measurements have been performed over 
a frequency range not limited to the adiabatic regime, so they can be compared with the 
prediction of the extended model proposed in Ref. [|]. 

The mirrors are part of two high-finesse Fabry-Perot interferometers whose length vari- 
ations are measured through phase- sensitive detection of resonant laser radiation. Our 
experimental apparatus was originally conceived to realize and test an optical readout for 



the gravitational wave bar detector AURIGA [11|; however we were able to observe the effect 



by actively modulating the intracavity light power, with minimal modifications of the setup. 
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FIG. 1: Scheme of the experimental setup: C#, Fabry-Perot cavity; 01, optical isolator; HW, half- 
wave plate; EOM^t, electro-optic modulator; P, polarizer; BS, beam-splitter; D#, photo-detector; 
M, mixer. Dashed lines mark out active thermal stabilization. 

The experimental setup is sketched in Fig. [|. The laser system and the detection are 
described in details in Ref. fl2 |. The light source is a cw Nd:YAG laser based on a non-planar 



ring oscillator geometry (Lightwave, mod. 126-1064-100), emitting 100 mW at 1064 nm. A 
non-resonant electro-optic modulator (EOM1) followed by a polarizer acts as an amplitude 
modulator for the transmitted laser radiation: the EOM1 has one crystal axis rotated by an 
angle of about 0.12 rad with respect to the input (vertical) polarization while the polarizer 
restores the vertical polarization. Applying a voltage to the EOM1 changes the ellipticity 
of the beam exiting EOM1 and thus the amplitude of light transmitted by the polarizer. 
By properly setting the temperature of the EOM1 we can reach and maintain a working 



point for where the variation of the transmitted light is proportional to the applied voltage. 
With no signal fed to the EOM1 the transmission after the polarizer is about 92%. This 
setup is part of a power stabilization loop described elsewhere [|1^] but in this work it is used 
merely for power modulation. A second, resonant electro-optic modulator (EOM2) phase 
modulates the laser output at 13.3 MHz. The two EOMs are enclosed in an aluminum 
box whose temperature is actively stabilized at ~ 35°C within 0.04°C A polarizer after 
EOM2, parallel to the previous one, removes any ellipticity that might originate in the 
phase modulation. The beam then crosses a first 50% beam splitter: the transmitted beam 
is detected by a photodetector (Dl) for power monitoring; the reflected beam is successively 
divided in two by a second 50% beam splitter. Each of the outgoing beams, of about 10 mW, 
passes an optical circulator formed by a polarizing beam-splitter and a quarter-of-wave plate 
and is then directed towards a Fabry-Perot cavity (CI and C2 in Fig. 0). The reflected beams 
are detected by photodiodes and the signals are demodulated and filtered, according to the 



FM sideband technique [13 



The two Fabry-Perot cavities are made by mirrors coated on fused silica substrates, which 
are optically contacted to a 20 cm long Zerodur spacer. Each cavity is placed on a double 
stage mechanical suspension which isolates from mechanical noise and it is enclosed in a 
vacuum chamber evacuated by an ion pump. The cavity temperature is actively stabilized 
within 1 mK. By varying the cavity temperature between 20°C and 100°C, it is possible 
to change the length of the cavity by roughly 0.5 /im, corresponding to one Free Spectral 
Range (FSR), and consequently to tune coarsely the optical resonances of the cavities. Fine 
adjustments of the frequency of CI is accomplished by a system of piezoelectric actuators 
(PZT) and invar rods which compresses the Zerodur spacer. Main figures of the two cavities 
are summarized in Table |. These cavities are designed as length reference for the AURIGA 
optical readout which is optimized for a frequency range around 1 kHz where it is charac- 
terized by a relative displacement noise of few 10~ 18 m/v'Hz. This noise increases as 1/f 
for lower frequencies 

The FM sideband signal from C2 is used as discriminator for frequency locking the laser. 
Once two resonance peaks of the two cavities have been coarsely superimposed, the laser 
frequency is locked to a resonance peak of C2 with a servo loop which has a unit gain 
frequency of 30 kHz and a gain of 120 dB around 1 kHz. CI is then finely tuned by 
acting on the PZT. The FM sideband signal from this cavity is used to monitor the relative 



5 



TABLE I: Parameters of the two cavities. We report the beam waist for each mirror, as usually 
defined in Gaussian optics, which differ from the beam radius ro, as defined in the text, by a factor 
V2. 



Cavity 


CI 


C2 


Length (m) 


0.2 


0.2 


FSR (MHz) 


750 


750 


Finesse 


16000 


36000 


Intracavity power (W) 


20 


240 


Mirror 


in out 


in out 


Radius of curvature (m) 


1 1 


oo 1 


Waist (mm) 


0.336 0.336 0.368 0.412 



frequency change and thus relative variations of the cavity lengths. 

To observe the photothermal effect we apply a sinusoidal voltage at frequency u m = u m /27i 
to EOM1. The power incident on both cavities is thus modulated around an average value 
Pq, according to 

P = P + P m smu m t, 

where P m is the modulation depth. Consequently the light power absorbed by the mirrors is 
modulated and a cavity length change is originated due to the non-zero thermal expansion 
coefficient of the fused silica substrates. We simultaneously acquire the signal from Dl, 
to monitor the power modulation, and the FM sideband signal from CI; these two signals 
are successively phase-detected at the frequency is m with a digital, numerical lock-in. The 
modulation frequency u m ranges between 5 x 10~ 3 Hz and 5 x 10 2 Hz, with constant modu- 
lation depth (P m /Po ~ 0.01). Measurements have also been made by fixing the modulation 
frequency at u m = 0.05 Hz and varying the modulation depth. 

If the two cavities were exactly identical we would observe no effects, since each cavity 
would experience the same length variation. However, as shown in Table |, due to different 
mirror losses the two cavities have different finesse and intracavity power, whereas the im- 
pinging laser power is the same. In particular, the power circulating into C2 is one order 
of magnitude larger than that into CI. Hence, we expect a differential length variation of 
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FIG. 2: Differential length variation 51 as a function of the relative modulation amplitude P m /Po, 
for a fixed modulation frequency v m = 0.05 Hz. The error bars are due to the uncertainty of the 
slope of the Pound-Drever signal, which calibrates the length scale. 

amplitude 81 between the two cavities and, as will be explained later, we can reasonably 
assume the observed effect is mainly due to C2. 

We have first checked that the observed effect is proportional to the modulation amplitude 
P m of power incident on the cavities. In Fig. |] the experimental differential length changes 
SI are plotted as a function of the relative modulation amplitude P m /Po, with v m = 0.05 Hz. 
The linear fit has a reduced \ 2 °f 0.75. 

In Fig. ^ we plot the experimental data for 51 as function of the modulation frequency 
u m . The reported error bars results from averaging over repeated acquisitions. These data 
have been fitted according to the model presented in Ref. 0. From Eqs. (|2]) and (|j) we 
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obtain: 

6l{u) = L oy jK{u/u c ), (4) 

where 

L = -(l + a)-SP ahs , (5) 

7T K 

and 5P a bs is the difference in the modulated power absorbed by the mirrors of the two 
cavities. The resulting fitting curve, also displayed in Fig. |3|, shows an excellent agreement 
between theory and experimental data (reduced x 2 = 0.5). /,From the fit we get the scale 
factor Lq = (3.3 ± 0.4) x 10~ 13 m and the cutoff frequency v c = 2.8 ± 0.6 Hz, where the 
reported error is the standard deviation. ^From Eq. ([!]) we calculate for v c = Lo c j2i\ the 
values of 2.7 Hz, 2.3 Hz and 1.8 Hz using respectively the three beam waists reported in 
Table | and the following figures for fused silica at room temperature ||: k = 1.4 W/m K, 
p = 2.2 x 10 3 kg/m 3 , C = 6.7 x 10 2 J/kg K fll5]| . We notice that the beam radius r , as 
defined above, is a factor l/y/2 with respect to the mirror beam waist as usually defined 
in Gaussian optics, which we report in Table | for each mirror. The agreement with the 
experiment is very good. To emphasize the peculiarity of the model we adopted, in Fig. [3] we 
also plot the best fitting curve (reduced \ 2 — 8) °f the experimental data with a single-pole, 



low-pass function, 51 = Lq/J 1 + {v/v c ) 2 . The functional form of 51 is a striking argument 
confirming the calculations of Cerdonio et al. 

As far as the size of the effect is concerned, we actually ignore the fractional amount of 
power which is absorbed at each reflection. From the reflection and transmission coefficients 
of each cavity we calculate the mirror losses, however these losses are not only due to 
absorption but also to scattering. Anyhow, we can give an estimate of the mirrors absorption 
coefficient, assuming that it is roughly the same for all the four mirrors, and that the losses 
are mostly due to scattering. From the parameter L resulting from the fit and Eq. (||), 
with a = 5.5 x 10~ 7 K _1 and a = 0.17 for fused silica, we get <5P a bs = 2.3 x 10~ 6 W for the 
modulated absorbed power. The differential intracavity power 5P C is 220 W, as from the 
value reported in Table |, and about 1% of it is actually modulated. Dividing the absorbed 
modulated power by the intracavity modulated power and by the number of mirrors (2) we 
obtain a mean value for the absorption coefficient of each mirror of about 5 x 10~ 7 , which 
is quite reasonable and in agreement with the estimate reported in Ref. 0. 

In conclusion we have experimentally investigated the photothermal effect on the effective 
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FIG. 3: Differential length variation 51 as m function of th.6 modulation frequency i^rw 

The solid 

line is the fitting curve of the experimental data with the expression of Eq. (|j). The dashed line is 
the best fit of the experimental data with a single-pole, low-pass function. 

length of high finesse Fabry-Perot interferometers. The measurements have been performed 
by modulating the power impinging on two cavities. The differential length change was 
measured for different modulation depths and frequencies, by means of an actively stabilized 
laser system with low residual frequency fluctuations and thanks to low noise, acoustically 
and thermally isolated cavities. The frequency functional dependence of the effect is per- 
fectly explained by the model proposed by Cerdonio et al. and the experimental parameters 
obtained from the fitting procedure are in very good agreement with independent estima- 
tions. Besides the interest for the physical phenomenon itself, this experimental study is a 
step toward a full understanding of thermal effects which are of critical importance for high 
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sensitivity displacement measurements. 
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